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ABSTRACT 

We have developed and successfully tested a new self-consistent method to reliably identify pre-main se- 
quence (PMS) objects actively undergoing mass accretion in a resolved stellar population, regardless of their 
age. The method does not require spectroscopy and combines broad-band V and / photometry with narrow- 
band Ha imaging to: (1) identify all stars with excess Ha emission; (2) convert the excess Ha magnitude into 
Ha luminosity L(Ha); (3) estimate the Ha emission equivalent width; (4) derive the accretion luminosity L acc 
from L(Ha); and finally (5) obtain the mass accretion rate M acc from L acc and the stellar parameters (mass and 
radius). By selecting stars with an accuracy of 15 % or better in the Ha photometry, the statistical uncertainty 
on the derived M acc is typically < 17% and is dictated by the precision of the Ha photometry. Systematic 
uncertainties, of up to a factor of 3 on the value of M acc , are caused by our incomplete understanding of the 
physics of the accretion process and affect all determinations of the mass accretion rate, including those based 
on a spectroscopic Ha line analysis. 

As an application of our method, we study the accretion process in a field of 9.16 arcmin 2 around SN 1987 A, 
using existing Hubble Space Telescope photometry. We identify as bona-fide PMS stars a total of 133 objects 
with a Ha excess above the 4a level and a median age of 13.5 Myr. Their median mass accretion rate of 
2.6 x 1O~ 8 M0 yr" 1 is in excellent agreement with previous determinations based on the {/-band excess of the 
stars in the same field, as well as with the value measured for G-type PMS stars in the Milky Way. The accretion 
luminosity of these PMS objects shows a strong dependence on their distance from a group of hot massive stars 
in the field and suggests that the ultraviolet radiation of the latter is rapidly eroding the circumstellar discs 
around PMS stars. 

Subject headings: accretion, accretion disks - stars: formation - stars: pre-main-sequence - Magellanic Clouds 
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1. INTRODUCTION 

In the current star formation paradigm, low-mass stars grow 
in mass over time through accretion of matter from a circum- 
stellar disc (e.g. Lynden-Bell & Pringle 1974; Appenzeller 
& Mundt 1989; Bertout 1989). This accretion process is be- 
lieved to occur via the magnetic field lines that connect the 
star to its disc and that act as funnels for the gas (e.g. Ktinigl 
1991; Shu et al. 1994). The strong excess emission observed 
in most of these pre-main sequence (PMS) stars is believed to 
originate through gravitational energy, released by infalling 
matter, that ionises and excites the gas. Thus, the excess lu- 
minosity of these objects (the classical T Tauri stars) can be 
used to measure the mass accretion rate. 

A reliable measurement of the rate of mass accretion onto 
PMS stars is of paramount importance for understanding the 
evolution of both the stars and their discs (see e.g. review by 
Calvet et al. 2000). Of particular interest is determining how 
the mass accretion rate changes with time as a star approaches 
its main sequence, whether it depends on the final mass of 
the forming star and whether it is affected by the chemical 
composition and density of its parent molecular cloud or by 
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the proximity of hot, massive stars. 

Ground-based spectroscopic studies of nearby young star- 
forming regions (e.g. in Taurus, Auriga, Ophiuchus, Orion) 
show that the mass accretion rate appears to decrease steadily 
with time, from ~ 10~ 8 M Q yr" 1 at ages of ~ 1 Myr to < 
1O _9 M yr" 1 at - 10 Myr (Muzerolle et al. 2000; Sicilia- 
Aguilar et al. 2005; 2006). While at face value this is in line 
with the expected evolution of viscous discs (Hartmann et al. 
1998), the scatter on the data is very large, exceeding 2 dex 
at any given age. This is at least in part explained by the re- 
cent finding that the mass accretion rate seems also to depend 
quite steeply on the mass of the forming star. Muzerolle et al. 
(2003; 2005), Natta et al. (2004; 2006), White & Hillenbrand 
(2004) and Sicilia-Aguilar et al. (2006) consistently report 
that the accretion rate M acc decreases with the stellar mass 
M as M acc oc M 2 , albeit with large uncertainties on the actual 
value of the power-law index. On the other hand, Clarke & 
Pringle (2006) warn that such a steep decline might be spuri- 
ous and strongly driven by detection/selection thresholds. 

While observational uncertainties and potential systematic 
errors may well contribute to the present uncertainty, the true 
limitation in this field of research currently comes from the 
paucity of available measurements. Indeed, all the results so 
far obtained are based on the mass accretion rates of a small 
number of stars (just over one hundred), all located in nearby 
Galactic star forming regions, covering a very limited range of 
ages and no appreciable range of metallicity (all clouds men- 
tioned above having essentially solar metallicity; e.g. Padgett 
1996). 

The origin of this limitation is to be found in the technique 
used to measure PMS mass accretion rates. Regardless as 
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to whether the latter are derived from the analysis of veiling 
in the photosperic absorption lines (usually at ultraviolet and 
optical wavelengths) or through a detailed study of the pro- 
file and intensity of emission lines (usually Ha, Pa/3 or Br7), 
these thechniques typically require medium to high resolution 
spectroscopy for each individual object. Even with modern 
multi-object spectrographs at the largest ground-based tele- 
scopes, this approach is not very efficient and cannot reach 
beyond nearby star forming regions or the Orion cluster, as 
crowding and sensitivity issues rule out but the brightest ob- 
jects in more distant regions like the centre of the Milky Way 
or the Magellanic Clouds. While in principle the PMS mass 
accretion rate can be estimated from the U -band excess if the 
spectral type is known (see e.g. Gullbring et al. 1998; Ro- 
maniello et al. 2004), this method is potentially subject to 
large uncertainties in the extinction and in the determination 
of the correct spectral type and effective temperature when 
spectroscopy is not available. 

In order to overcome these limitations and significantly ex- 
tend the sample of stars with a measured mass accretion rate, 
we have developed and successfully tested a new method to 
reliably measure M acc that does not require spectroscopy. The 
method, described in the following sections, combines broad- 
band V and / photometry with narrow-band Ha imaging and 
allows us to identify all stars with excess Ha emission and to 
derive from it the accretion luminosity L acc and hence M acc for 
hundreds of objects simultaneously. 

The paper is organised as follows: in Section 2 we ad- 
dress the identification of PMS stars via their colour excess 
and in Section 3 we show how to reliably convert the mea- 
sured excess into Ha luminosity and equivalent width. Sec- 
tion 4 shows how the accretion luminosity and mass accre- 
tion rate can be obtained from the Ha luminosity and Sec- 
tions presents an application of the method to the field of 
SN 1987A, in the Large Magellanic Cloud (LMC). A thor- 
ough discussion of all systematic uncertainties involved in the 
method is provided in Section 6, while a general discussion 
and conclusions follow in Section 7. The Appendix provides 
additional figures and tabular material. 

2. IDENTIFICATION OF PMS STARS THROUGH THEIR COLOUR 
EXCESS 

One of the characteristic signatures of the accretion pro- 
cess on to PMS stars is the presence of excess emission in 
Ha (see Calvet et al. 2000). Although surveys to search for 
stars with Ha excess emission over large sky areas have been 
traditionally carried out via slitless spectroscopy (e.g. Her- 
big 1957; Kohoutek & Wehmeyer 1999), extensive studies 
in which narrow-band Ha imaging is combined with broad- 
band photometry also exist (e.g. Parker et al. 2005; Drew et 
al. 2005). 

Recently, Romaniello et al. (1998), Panagia et al. (2000) 
and Romaniello et al. (2006) used HST/WFPC2 photometry 
in Ha (F656N) and the R band (F675W) to provide the first 
direct detection of extra-galactic PMS stars with strong Ha 
emission in the LMC. Their approach is based on the widely 
adopted use of the R band magnitude as an indicator of the 
level of the photospheric continuum near the Ha line. There- 
fore, stars with strong Ha emission will have a large R-Ha 
colour. From the analysis of the specific F656N and F675W 
filter response curves of the WFPC2 instrument, Panagia et 
al. (2000) conclude that an equivalent width (W eq ) of the Ha 
emission line in excess of 8 A will result in a colour excess 
R-Ha> 0.3. 



This way of identifying PMS stars is more accurate and reli- 
able than the simple classification based on the position of the 
objects in the Hertzsprung-Russell diagram, i.e. stars placed 
well above the main sequence (e.g. Gilmozzi et al. 1994; 
Hunter et al. 1995; Nota et al. 2006; Gouliermis et al. 2007). 
The reason is twofold: ( i) contamination by older field stars 
and the effects of differential extinction may lead to an over- 
estimate of the actual number of candidate young, low-mass 
PMS stars, i.e. those "above" the MS; ( ii) older PMS stars, al- 
ready close to their MS, cannot be properly identified with this 
method, thereby resulting in an underestimate of the number 
of stars that are still forming. When PMS stars are identified 
by means of their R-Ha colour excess, both problems disap- 
pear. 

Similarly, this approach is more practical, efficient and eco- 
nomical than slitless spectroscopy for studying PMS stars 
in dense star forming regions in the Milky Way and nearby 
galaxies, which are today accessible thanks to high-resolution 
observations, in particular those made with the HST. 

We will show in Section 3 that the level of the photospheric 
continuum near Ha can be directly derived from the observed 
V and / magnitudes. This means that PMS stars with an 
Ha excess can be identified without actually measuring the 
canonical R- Ha index. On the other hand, since R-Ha has 
traditionally been used in PMS stars studies, we first show 
here how one can derive analytically the R magnitude of an 
object from the measured V and I magnitudes. In this way, 
we intend to make it possible for those who do not dispose 
of R-band photometry for their fields to derive it in an accu- 
rate way by using measurements in the neighbouring V and I 
bands. The latter are more commonly used than the R band 
in photometric studies of resolved stellar populations and are 
particularly popular among HST observers. 

2.1. Deriving Rfrom V and I 

The absence of remarkable spectral features between 
5 000 A and 9 000 A in the stellar continuum of dwarfs and 
giants with effective temperatures in the range 4000K < 
7k < 10000K implies that their V-R and V -I colours de- 
pend primarily on the effective temperature of the stars and 
less on their metallicity. In effect, V-I remains a useful 
colour index for temperature determinations (von Braun et al. 
1998; Bessell, Castelli & Plez 1998). 

The colour-colour sequences in the UBVRIJKLMN pass- 
bands published by Johnson (1966) revealed a very tight cor- 
relation between the mean V-R and V-I colours of local 
dwarfs and giants (see Figure[T|i. The correlation is remark- 
ably linear in the range < V-I < 2 (i.e. 4000K < T efi < 
10000 K), where (V-R) = 0.58 x (V-I) (dashed line in Fig- 
ure[TJ, and it remains linear for 2 < V-I < 4 (i.e. 3 000K 
< 7k < 4000K), albeit with a shallower slope (0.45; dot- 
dashed line). It is, therefore, possible and meaningful to at- 
tempt to derive the value of the R magnitude from the knowl- 
edge of V and /, particularly because no extrapolation beyond 
the available wavelength range is involved. The simple inver- 
sion of the relationship above gives R = 0.42 V + 0.58/ for the 
system of Johnson0 

Interestingly, this relationship is very similar to the one 
that we would obtain if we assumed that the proportion with 

6 We note here that Johnson's / is not the same as the much more com- 
monly used Cousin's /. Therefore, the colour relationships given here cannot 
be used for the Cousin system. Proper relationships for the most common 
magnitude systems are provided in the Appendix. 
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FIG. 1. — The tight correlation between V — R and V-I from the colour- 
colour sequences of Johnson (1966). Circles indicate dwarfs and crosses gi- 
ants. 

which the V and / magnitudes contribute to the R magnitude 
is simply dictated by the ratio of the effective wavelengths 
of the filters. The effective wavelengths of the V, R and I 
bands in Johnson's system are, respectively, \ v e = 5500 A, 
X e R = 7000 A and = 9000 A , which imply that the V band 
should contribute for (7000-5500)/(9000-5500) ~ 0.43 and 
the I band for (9000 - 7000)/(9000 - 5500) ~ 0.57 to the R- 
band magnitude. These values are remarkably close to those 
observed. 

The apparent disarming simplicity of this conclusion, how- 
ever, should not lead one to think that this is the case for any 
photometric system. In particular, for spectral bands much 
wider than Johnson's, such as some of the HST/WFPC2 and 
HST/ACS filters, the matter becomes more complicated be- 
cause the true effective wavelength can depend significantly 
on the spectral properties of the source (i.e. its temperature). 
In these cases, the relationship should be derived via direct 
observations in the specific bands or through synthetic pho- 
tometry, when the properties of the photometric system are 
well caracterised. 

Due to the large number of inter-related observing modes 
offered by the scientific instruments on board the HST, syn- 
thetic photometry has long been established as the most prac- 
tical and efficient approach to their calibration (see Koornneef 
et al. 1986; Home 1988; Sirianni et al. 2005). For this reason, 
reference tables exist that accurately describe the main optical 
components of the telescope and instruments as well as the re- 
sponse and quantum efficiency of the detectors as a function 
of wavelength. Therefore, thanks also to the availability of 
standardised software packages (e.g. Synphot; see Laidler et 
al. 2008), it is relatively easy to derive synthetic magnitudes 
from observed or theoretical stellar spectra for various HST 
instrumental configurations. 

We have used both observed and theoretical stellar spectra, 
taken respectively from the Bruzual-Persson-Gunn-Stryker 
(BPGS) Spectrophotometry Atlas (see Gunn & Stryker 1983) 
and from the stellar atmosphere models of Bessell et al. 
(1998), to compute colour relationships between R-like mag- 
nitudes and V-/-like colours for various HST instruments and 
ground photometric systems. Details on how the colour trans- 



FlG. 2. — Colour-colour diagram of all stars in the SN 1987A field with 
combined mean photometric error £3 < 0.05 mag. The dashed line shows 
the colour relationship derived for the bands in question (F555W, F675W 
and F814W) from the models of Bessell et al. (1998) and is in excellent 
agreement with the observations. 

formations were derived and a table of coefficients (TableO 
can be found in the Appendix, while a brief comparison with 
HST data is given in the following section. 

2.2. Comparison with the data 

To verify the validity of the predicted colour relationships, 
we have used WFPC2 photometry of a field around SN 1987 A 
studied by Romaniello (1998), Panagia et al. (2000) and Ro- 
maniello et al. (2002). The observations cover a field of 
9.16arcmin 2 and were collected over three epochs, namely 
1994 September, 1996 February and 1997 July. The advan- 
tage of this catalogue is that the magnitude of each star in it 
is individually corrected for the effects of interstellar extinc- 
tion by using the information provided simultaneously by all 
broad-band filters (see Romaniello et al. 2002 for more de- 
tails). While it would be possible to include the effects of 
reddening in the calculations of the colour relationships of 
Table[3] it is more appropriate in this case to correct the ob- 
served magnitudes, since considerable differential extinction 
is likely to be present in this field (see Panagia et al. 2000; De 
Marchi & Panagia, in preparation). 

From the catalogue of Romaniello et al. (2002), we have 
selected all those stars whose mean error S3 in the three bands 
F555 W, F675 W and F8 14W does not exceed 0.05 mag, where 



fr = A±fc& (1) 

and (5 5 55, <5 67 5 and 5%\n are the photometric uncertainties in 
each individual band after correction for interstellar extinc- 
tion. A total of 9 635 stars satisfy the condition set by Equa- 
tionQ] out of 21 955 objects in the complete catalogue. 

The colour-colour diagram for the stars selected in this 
way is shown in Figure[2] The dashed line in that figure 
corresponds to the colour relationship derived for the bands 
in question from the models of Bessell et al. (1998; see 
the Appendix) and is practically indistinguishable from the 
best linear regression fit to the data, which has a slope of 
0.645 ± 0.003. Although these models were selected to match 
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FIG. 3. — Comparison between the measured and derived R — Ha colour excess. While R denotes the magnitude actually measured in the F675W band, R* 
indicates the magnitude derived from V (F555W) and / (F814W) by using the relationship of Table[5] Diamonds represent the measured average R-Ha colour 
of the population while the dashed lines is the colour predicted my theoretical models. Thick dots corresponds to stars with R—Ha excess larger than 4cr. See 
text for an explanation of the remaining symbols. 



the metallicity of the LMC, having used the relationship de- 
rived with the BPGS atlas or that for [M/H] = -1 would have 
still resulted in a very good fit, within the observational un- 
certainties. This result confirms that, at least over the range 
explored here, metallicity does not play a dominant role in 
the colour relationships. 

We provide a further proof of the reliability and accuracy 
of the interpolated R-band magnitudes in Figure[3] where we 
compare the identification of PMS stars based on the mea- 
sured and interpolated R - Ha colour. To distinguish the ob- 
served from interpolated R-band magnitude, we indicate the 
latter with R*, as explained belowQ In this case, from the 
photometry of Romaniello et al. (2002) already corrected for 
reddening, we have further restricted the selection to objects 
with S3 < 0.02, comprising 3241 stars. 

Diamonds in Figure[3] represent the average R-Ha colour 
of the population and, as such, define the reference with re- 

7 Hereafter, we will use for simplicity the symbols V, R, I and Ha to in- 
dicate the respective WFPC2 bands. When differences between photometric 
systems are important (e.g. in TablefSJ, all bands will be indicated with their 
specific name. 



spect to which one should look for excess emission. The 
dashed lines represent the theoretical colour relationship ob- 
tained using the Bessell et al. (1998) model atmospheres 
mentioned above for normal main-sequence stars and is in 
rather good agreement with the average observed colours in 
the range -0.1 < V— / < 1 (filled diamonds). The rms devia- 
tion between the model and the data over this range amounts 
to 0.03 mag and is dominated by the systematic departure in 
the range 0.2 < V— / < 0.7. The latter could be due to a num- 
ber of reasons, including the inability of the models to accu- 
rately describe the properties of the Ha line since they are 
meant to be used for broad-band photometry. The departure 
of the models from the data at V-l > 1 stems from both small 
number statistics and the fact that the bulk of the population 
there comes from red giants, which have lower surface grav- 
ity and possibly lower metallicity than those assumed in the 
models. In the following, we will concentrate our analysis to 
the range -0.1 < V — / < 1 marked by the filled diamonds. 

Figure^ is based on the observed data and we indicate as 
large dots all stars whose observed R-Ha colour exceeds 
the local average by at least four times the individual pho- 
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tometric uncertainty in the colour, i.e. the one specific to each 
star. Consequently, these objects must be regarded as hav- 
ing a bona-fide Ha excess above the 4 a level. There are 199 
such objects in panel (a). In panel (b) we have replaced the 
observed R magnitude with its interpolated value R*. All but 
16 of the 199 stars with 4 a excess in panel (a) have a 4 a 
excess also in panel (b) and are indicated as large dots. The 
remaining 16 objects, marked as triangles, fail this require- 
ment as their colour is closer to the local average colour than 
four times their specific combined photometric uncertainty. 

In panels (c) and (d) we study the uncertainty in the other di- 
rection, i.e. from the interpolated to the observed plane. Panel 
(c) simulates the case in which, no /?-band data being avail- 
able, one uses the interpolated R* value from Table[3]to iden- 
tify stars with Ha excess emission. Stars indicated as large 
dots in panel (c) are those with interpolated R* -Ha colour 
exceeding the local average by at least four times the indi- 
vidual photometric uncertainty in the colour, i.e. those with 
Ha excess above the 4cr level. A total of 194 stars meet this 
condition. All but 1 1 of them would have met the 4 a excess 
condition if we had used the actually observed /?-band magni- 
tude. They are indicated as tlarge dots in panel (d), while the 
1 1 objects failing the condition are marked as triangles. 

From Figure|3jwe can draw some quantitative conclusions 
on the statistical significance of the detection of the colour ex- 
cess in Ha. First of all, it does not surprise that some objects 
detected as bona-fide Ha excess stars in panel (a) are rejected 
in panel (b), since the typical photometric uncertainty on R* 
is a factor of y/2 higher than that on R, as one must com- 
bine the uncertainties on V and / that concur to the estimate 
of R* (although the major contribution to the uncertainty on 
the R-Ha colour of a star comes from the Ha band). For the 
same reason, it is to be expected that panel (c) will yield less 
objects with bona-fide Ha excess than panel (a). 

Considering the number of detections and mismatches of 
Figure[3] we estimate that using the interpolated R* in place 
of R a fraction of ( 1 99 - 1 94)/ 1 99 = 2 .5 % of the objects with 
Ha excess at the 4 a level would be missing and that for about 
11/194 = 5.7% of the detections the assigned statistical sig- 
nificance would be higher than indicated by direct R measure- 
ments. We conclude that our interpolation scheme is quite re- 
liable in that it can provide the correct identification of PMS 
stars in no less 94 % of the cases. 

3. QUANTIFYING THE EXCESS Ha EMISSION: FROM COLOUR 
EXCESS TO LINE LUMINOSITY AND EQUIVALENT WIDTH 

While helpful to accurately identify PMS stars, the R—Ha 
colour excess alone does not immediately provide an absolute 
measure of the Ha luminosity L(Ha) nor of the equivalent 
width W eq (Ha), which are necessary to arrive at the mass ac- 
cretion rate. In order to measure L{Ha), one needs a solid es- 
timate of the stellar spectrum in the Ha band (i.e. without the 
contribution of the emission). Similarly, to measure W eq (Ha), 
a solid estimate of the level of the stellar continuum inside 
the Ha band is needed. Since the R band is over an order of 
magnitude wider than the Ha filter, it cannot provide this in- 
formation accurately, but only gives a rough estimate of the 
continuum level. In this section we present a simple method 
that allows us to measure L(Ha) for the objects with an Ha 
excess and to derive the level of their spectral continuum in- 
side the Ha band, determining in this way also W eq {Ha). 

3.1. The Ha line luminosity 



The method is based on the simple consideration that, at 
any time, the largest majority of stars of a given effective 
temperature T e ff in a stellar population will have no excess 
Ha emission. Therefore, for stars of that effective tempera- 
ture, the average value of a colour index involving Ha (for 
instance V -Ha, and thus not just R-Ha) effectively defines 
a spectral reference template with respect to which the Ha 
colour excess should be sought. Note that this is not only true 
for populations comprising both young and old stars, such as 
those typical of extragalactic stellar fields, but also for very 
young populations since PMS objects show large variations 
in their Ha emission over hours or days (e.g. Fernandez et al. 
1995; Smith et al. 1999; Alencar et al. 2001), with only about 
one third of them at any given time being active Ha emitters 
above a W eq (Ha) ~ 10 A threshold (Panagia et al. 2000). 

To better clarify the working of our method, we display in 
Figure|4]the V -Ha colour as a function of V-I for a set of 
stars taken from the catalogue of Romaniello et al. (2002). 
This time, we have selected a total of 4 156 stars, namely all 
objects with S3 < 0.1 mag, where the mean uncertainty £3 is 
defined as in Equation[T| but for the V, I and Ha bands instead 
of V, I and R. Therefore, no R-band information is used in 
this case and, in fact, no R-band data is needed. It should 
be noted that £3 is dominated by the uncertainty on the Ha 
magnitude, while the median value of the uncertainty in the 
other two bands is 8y = 0.014 and 5j = 0.016, respectively. 

The thick dashed line in FigureH] represents the average 
V - Ha colour obtained as the running median with a box- 
car size of 100 points. The thin solid line shows the colours in 
these filters for the model atmospheres of Bessell et al. (1998) 
mentioned above (note again the < 0.1 mag discrepancy be- 
tween models and observations around V-/ ~ 0.6, most likely 
due to the coarse spectral sampling of the Ha line in the mod- 
els). At the extremes of the distribution (V -I < -0.1 and 
V-I > 1) the density of observed points decreases consider- 
ably. To determine the reference template there, we have re- 
duced the box-car size to 10 points and averaged the running 
median obtained in this way with the theoretical models (thin 
solid line). The thick dot-dashed line in Figure|4] shows the 
best linear fit to the resulting average (which is extrapolated 
forV-/>2). 

A total of 189 objects, indicated by large dots in FigureH] 
have a V -Ha index exceeding that of the reference template 
at the same V-I colour by more than four times the uncer- 
tainty on their V -Ha values. These are the objects with a 
V —Ha excess at the 4 a level and, as we shall see later, most 
of them are bona-fide PMS stars. 

Since the contribution of the Ha line to the V magnitude is 
completely negligible, the magnitude AHa corresponding to 
the excess emission is simply: 

AHa = (V-Ha) ohs -(V-Ha) le{ (2) 

where the superscript obs refers to the observations and ref 
to the reference template. Once AHa is determined in this 
way, the Ha emission line luminosity L(Ha) can be imme- 
diately obtained from the photometric zero point and abso- 
lute sensitivity of the instrumental set-up and from the dis- 
tance to the sources. We have assumed a distance to the 
LMC and, more specifically, to SN 1987A of 51.4 ± 1.2kpc 
(Panagia et al. 1991, later updated in Panagia 1999), whereas 
the photometric properties of the instrument were taken from 
the WFPC2 Instrument Handbook (Heyer & Biretta 2004; 
namely inverse sensitivity PHOTFLAM= 1.461 x 10" 16 and 
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FIG. 4. — Colour-colour diagram of the selected 4 156 stars in the field of 
SN 1987A. The dashed line represents the running median V — Ha colour, 
obtained with a box-car size of 100 points, whereas the thin solid line shows 
the model atmospheres of Bessell et al. (1998). A total of 189 objects with a 
V —Ha excess larger than 4<r are indicated with large dots. 
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FIG. 5. — Histogram of the Ha luminosities of the 189 stars with Ha excess 
at the 4<r level (dashed line) and of for the 133 bona-fide PMS that also have 
an emission W C q(//a) > 10 A and are not Ae/Be stars (solid line). 

photometric zero-point ZP(Vegamag)= 17.564). We derive a 
median value of the luminosity of the 189 objects with Ha 
excess of 4 x 10 31 erg/s or 1O~ 2 L0. A histogram of the indi- 
vidual L(Ha) values is shown in Figure[5] (dashed line). 

It should be noted that, because of the width of the Ha fil- 
ter, AHa includes a small contribution due to [Nil] emission 
features at 6 548 A and 6 584 A. Generally, their intensities do 
not exceed 1 .2 % and 3.5 %, respectively, of the Ha line inten- 
sity, and in many instances the [Nil] doublet lines are much 
weaker, i.e. < 0.1 % of the Ha line intensity (e.g. Edwards et 
al. 1987; Hartigan, Edwards & Ghandour 1995). These inten- 
sities should be even lower for PMS stars in the LMC because 
of the lower metallicity relative to the Milky Way. However, 
to take a conservative approach, we have simply adopted an 



average value that is half of the maximum for Galactic T-Tauri 
stars and we give as uncertainty half of the measured range. 
We will therefore adopt average values of (1.8 ± 1.8)% and 
(0.6 ± 0.6) % relative to the Ha intensity, for the 6 584 A and 
6 548 A [Nil] lines, respectively (note that the uncertainties 
after the ± sign here are meant to show the entire span of 
the range, not just the 1 a value). Considering the width of 
the Ha filter of the WFPC2, only the 6 548 A [Nil] line con- 
taminates the Ha measurements, whereas for measurements 
made with the Ha filter on board the ACS, both lines will be 
included (albeit only in part for the 6 548 A line). Although 
the resulting corrections on the Ha intensity are quite small 
(namely 0.994 for the WFPC2 F656N filter and 0.979 for the 
ACS F658N filter, on average), they are systematic and it is 
therefore good practice to apply them in all cases. 

Typically, the combined total uncertainty on our L(Ha) 
measurements is ~ 15 % and is completely dominated by the 
statistical uncertainty on the Ha photometry. The systematic 
uncertainty on the distance to the field of SN 1987 A accounts 
for 5 % (see Panagia 1999), whereas the typical uncertainty 
on the absolute sensitivity of the instrumental setup is of or- 
der 3%. 

3.2. The Ha equivalent width 

If the stars defining the reference template had no Ha ab- 
sorption features (we have established above that they do not 
have Ha in emission), their V - Ha index would correspond 
to that of the pure continuum and could then be used to derive 
the equivalent width of the Ha emission line, W eq (Ha). While 
this approximation might be valid for stars with a conspicu- 
ous Ha emission, it is clearly not applicable in general. It is, 
however, possible in all cases to derive the level of the Ha 
continuum by using properly validated models atmospheres. 

As discussed above, the models of Bessell et al. (1998) re- 
produce quite reliably the observed broad-band colours (see 
e.g. Figure|2]l. This means that, even though the mod- 
els might lack the resolution necessary to realistically repro- 
duce spectral lines (hence the small discrepancy between the 
dashed line and the squares in Figure[3]), the level of the con- 
tinuum can be trusted. One can, therefore, fit the contin- 
uum in those models in a range containing the Ha line (e.g. 
6 500 -6 620 A) and fold the resulting "Ha-less" model spec- 
tra through the instrumental set-up (e.g. with Synphot), so as 
to derive the magnitude of the sole continuum in the Ha band 
(Ha c ) as a function of spectral type or effective temperature. 
We present in the Appendix the relationships between V, I 
and Ha c for several HST instruments and filters. The differ- 
ence between the observed Ha magnitude and Ha € provides 
a direct measure of W eq (Ha), as we explain below. 

We recall that the equivalent width of a line is defined as: 



W eq = / (l-Px)d\ 



(3) 



where Pa is the profile of the line, or the spectrum of the 
source divided by the intensity of the continuum, and the 
integration is extended over the spectral region correspond- 
ing to the line. If the line profile is very narrow when com- 
pared to the width of the filter, i.e. if the line falls completely 
within the filter bandpass as is typical of emission lines, then 
W eq (Ha) is simply given by the relationship: 



W eq (Ha) = RWx [1-10 



-0.4x(ffQ-ffa')i 



(4) 
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FIG. 6. — Equivalent width of the selected 4156 stars in the field of 
SN 1987A, as a function of their V-l colour. Of the 189 stars with Ha excess 
at the 4<t level (circles and triangles), a total of 164 have W eq (//a) > 10A. 
While objects bluer than V— I ~ 0.2 are probably Ae/Be stars, there are 133 
redder stars with significant Ha excess (marked as circles) that we consider 
bona-fide low-mass PMS stars. 

where RW is the rectangular width of the filter (similar in 
definition to the equivalent width of a line), which depends 
on the characteristics of the filter. Values of RW for the Ha 
bands of past and current HST instruments are given in the 
Appendix (Table|4j>. As in the case of L(Ha), the statistical 
uncertainty on the equivalent width is dominated by the un- 
certainty on the Ha magnitude, typically of order 15 % in our 
selection. Moreover, the value of W eq {Ha) obtained in this 
way is also subject to some systematic uncertainties produced 
by the model atmospheres used to interpolate the level of the 
continuum, such as e.g. metallicity, surface gravity or spec- 
tral resolution, and should therefore be taken with caution. 
On the other hand, we show in the Appendix that the equiva- 
lent widths (in absorption) obtained by applying this method 
to the high-resolution model spectra of Munari et al. (2005) 
are in excellent agreement with those obtained via standard 
spectro-photometry (see Figure[T8l). 

Similarly, the values of W eq (Ha) that we obtain in the field 
of SN 1987 A are fully consistent with those of young PMS 
stars. In Figure[6] we show the value of W eq (Ha) as a func- 
tion of V— I for the 4 156 selected stars. As mentioned earlier, 
there are 189 objects with a AHa excess above the 4er level 
and they are indicated as thick symbols (circles and triangles) 
in Figure[6] Since the equivalent width shown in Figure|6] is 
that of the pure emission component, the spectra of stars with 
small Weq, say < 10 A, actually show an Ha in absorption. 
For this reason, we conservatively ignore stars with measured 
Weq(Ha) < 10 A, since this value is about the largest absorp- 
tion equivalent width expected for normal stars (see FigurefhSl 
in the Appendix). A total of 164 objects satisfy this condition, 
but, in light of our conservative threshold, this must be con- 
sidered as a lower-limit to the number of objects with genuine 
Ha excess. 

Values of W eq (Ha) for the sample range from 10A to 
650 A, with a median of 19 A. These values are in excellent 
agreement with those obtained by Panagia et al. (2000) from 
the R-Ha colour excess of the stars in this field and represent 
a refinement of those measurements, since these are based 



on a more accurate determination of the level of the contin- 
uum. These equivalent widths are typical of PMS stars (e.g. 
Muzerolle, Hartmann & Calvet 1998) and about two orders of 
magnitude smaller than those expected for pure nebular emis- 
sion (see discussion in Section 6). While objects bluer than 
V -I ~ 0.2 are probably Ae/Be stars, there are 133 redder 
stars with significant Ha excess (marked as circles in Fig- 
ure[6]) that are most likely bona-fide low-mass PMS stars of 
various ages, as we will show in the following sections. A 
histogram of their Ha luminosities is shown in Figure|5](solid 
line). 

4. ACCRETION LUMINOSITY AND MASS ACCRETION RATE 

It is generally accepted (see e.g. Hartmann et al. 1998) 
that the energy L acc released by the magnetospheric accretion 
process goes towards ionising and heating the circumstellar 
gas. In this hypothesis, the bolometric accretion luminosity 
L acc can be determined from a measurement of the reradiated 
energy and, in particular, via the measured Ha luminosity that 
is produced in the process. 

In order to determine the exact relationship between L(Ha) 
and L acc , we have used literature values of L acc and L(Ha) 
measurements as recently summarized by Dahm (2008) for 
PMS stars in the Taurus-Auriga association. As pointed out 
by Dahm himself, there may be a serious problem in de- 
termining an empirical relationship between L acc and L(Ha) 
measurements because the two quantities were determined 
from non-simultaneous observations. The intensity of the Ha 
line from PMS sequence stars is known to vary by about 20 % 
on a time scale of a few hours and by as much as a factor of 
2 - 3 in a few days (e.g. Fernandez et al. 1995; Smith et 
al. 1999; Alencar et al. 2001). For this reason, a relation- 
ship between L(Ha) and L acc is necessarily rather uncertain, 
although a very clear trend can be seen in the L acc vs. L(Ha) 
plot shown by Dahm (2008). His logarithmic best fit would 
provide a slope of 1.18 ±0.26 for such a relationship. On 
the other hand, theoretical models (e.g. Muzerolle, Calvet & 
Hartmann 1998) would predict logarithmic slopes of about 
unity for low accretion rates, hence for faint Ha luminosities, 
and shallower slopes for higher luminosities. 

In the absence of compelling evidence in favour or against 
a steep slope, we will adopt a logarithmic slope of unity for 
the empirical L acc vs. L(Ha) relationship, i.e. a constant ratio 
Lncc/L(Ha), and determine the proportionality constant from 
an elementary fit of the data summarised by Dahm (2008). On 
this basis we obtain: 



logL acc = ( 1 .72 ± 0.47) + \ogL(Ha) 



(5) 



As anticipated, the resulting uncertainty on L acc is rather 
large, as high as a factor of 3, but this is the best one can 
obtain with the available data sets. In view of the conve- 
nience and the feasibility of determining L acc from L(Ha) for 
PMS outside our own Galaxy, it would be of paramount im- 
portance to obtain a proper calibration of the L acc vs. L(Ha) 
relationshipby making simultaneous spectral observations of 
PMS stars over a suitably wide wavelength baseline, i.e. cov- 
ering at least the range 3000-8000 A, and possibly extending 
it to the near IR with measurements of Paschen and Brack- 
ett lines. High quality data of this type are already available 
for the Orion Nebula Cluster, collected by dedicated projects 
both from ground based observatories and with the HST (i.e. 
the Legacy Program on Orion; P.I. Massimo Robberto) and 
their analysis is in progress to provide a more accurate cal- 
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ibration of the L acc vs. L(Ha) relationship (Da Rio et al., in 
preparation; Robberto et al., in preparation; Panagia et al., in 
preparation). 

The median value of the accretion luminosity thus obtained 
for our sample of 133 PMS stars is ~ O.57L , in good agree- 
ment with the value of 0.54 L Q derived from the U-band ex- 
cess measured by Romaniello et al. (2004) for stars in the 
same field, especially considering that their selection of stars 
with excess emission is less stringent than our 4 a level (we 
discuss this point in Section 7). 

Once L acc is known, the mass accretion rate follows from 
the free-fall equation that links the luminosity released in the 
impact of the accretion flow with the rate of mass accretion 
M acc according to the relationship: 



GM* M ac 
R* 



(6) 



where G is the gravitational constant, and R* the mass 
and photospheric radius of the star and Ri n the inner radius of 
the accretion disc. The value of R in is rather uncertain and 
depends on how exactly the accretion disc is coupled with the 
magnetic field of the star. Following Gullbring et al. (1998), 
we adopt R m = 5R* for all PMS objects. 

As regards we derive it from the luminosity and ef- 
fective temperature of the stars, which in turn come from 
the observed colour (V-I) and magnitude (V), properly cor- 
rected for interstellar extinction as provided by Romaniello 
(1998), Panagia et al. (2000) and Romaniello et al. (2002). 
The adopted distance to SN 1987 A is 51.4kpc (Panagia et 
al. 1991; Panagia 1999), as mentioned above. The stel- 
lar mass M* was estimated by comparing the location of a 
star in the Hertzsprung-Russell (H-R) diagram with the PMS 
evolutionary tracks. As for the latter, we adopted those of 
DeglTnnocenti et al. (2008; see also Tognelli, Prada Moroni 
& DeglTnnocenti, in preparation), for metallicity Z = 0.007 or 
about one third Z Q , as appropriate for the LMC (e.g. Hill, An- 
drievsky, & Spite 1995; Geha et al. 1998). These new PMS 
tracks were specifically computed with an updated version of 
the FRANEC evolutionary code (see Chieffi & Straniero 1989 
and DeglTnnocenti et al. 2008 for details and Cignoni et al. 
2009 for an application). However, in order to assess how 
differences in the evolutionary models affect our results, in 
Section 6.1 we also consider PMS tracks from other authors, 
namely those of D'Antona & Mazzitelli (1997) and Siess, Du- 
four & Forestini (2000). 

Combining Equations[5]and[6] we can link the mass accre- 
tion rate M acc , in units of M© yr" 1 , to L(Ha): 



log 



Moyr" 1 



-7.39 + log ^ + log ^- log ^ 

S L© & R G S M© 



(7) 



L(Ha) R* AL 
= (-5 .67 ± 0.47) + log + log -1 - log -1 

Lq Kq Mq 



The distribution of mass accretion rates found in this way 
for the 133 PMS stars in our sample is shown in Figure|7] The 
median value of 2.9 x 10~ 8 M© yr" 1 is in remarkably good 
agreement with the median rate of 2.6 x 10~ 8 M© yr" 1 found 
by Romaniello et al. (2004) from the U-band excess. 

The uncertainty on M acc is dominated by the approximate 
knowledge of the ratio of L acc and L(Ha) in Equation[5] 
As explained above, empirical measurements and theoretical 
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FIG. 7. — Distribution of the mass accretion rates of the 133 bona-fide PMS 
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models suggest a value of L acc /L{Ha) ~ 52 but with an un- 
certainty of a factor of three. Since the ratio, however, is the 
same for all stars, the comparison between different objects 
is not hampered by this uncertainty, as long as the statistical 
errors are small. 

As for the other quantities in Equation 7, we discuss here 
the sources of statistical uncertainty, while systematic errors 
are addressed separately in Section 6. With our selection cri- 
teria, the typical uncertainty on L(Ha) is 15% and is dom- 
inated by random errors. The uncertainty on R* is typically 
7 %, including a 5 % systematic uncertainty on the distance 
modulus. As for the mass M*, since it is determined by com- 
paring the location of a star in the H-R diagram with evolu- 
tionary tracks, both systematic and statistical uncertainties are 
important. The uncertainty on the temperature is mostly sta- 
tistical and of order 3 %, while that on the luminosity is 7%, 
comprising both random errors (1 % uncertainty on the bolo- 
metric correction and 3 % on the photometry) and systematic 
effects (5 % uncertainty on the distance modulus). When we 
interpolate through the PMS evolutionary tracks to estimate 
the mass, the uncertainties on T e g and L imply an error of 
~ 7% on M*. However, an even larger source of systematic 
uncertainty on comes from the evolutionary tracks, as we 
explain in Section 6. In summary, the combined statistical un- 
certainty on M acc is 17 %. 

5. AGES OF PMS STARS 

Having identified a population of bona-fide PMS stars, with 
well defined accretion luminosities and mass accretion rates, 
it is instructive to place them in the H-R diagram. We do so in 
Figure[8]where bona-fide PMS stars are marked as large dots, 
with a size proportional to their M acc value, according to the 
legend. Also shown in the figure are the theoretical isochrones 
from the FRANEC models of DeglTnnocenti et al. (2008) for 
Z = 0.007 and ages as indicated (in units of Myr, from the 
stellar birth line; see Palla & Stahler 1993). The dashed line, 
corresponding to a 40 Myr isochrone, defines in practice the 
zero-age main sequence (MS). 

It is noteworthy that the majority of PMS objects in Fig- 
ure[8] are rather close to the MS and would have easily been 
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FIG. 8. — Location of the bona-fide PMS candidates in the H-R diagram, 
compared with the isochrones of Degl'Innocenti et al. (2008) for metallicity 
Z = 0.007 (the ages, in Myr, are indicated next to each isochrone). The size 
of the symbols is proportional to the mass accretion rate M acc (see legend). 

missed if no information on their Ha excess had been avail- 
able. Indeed, it is customary to identify PMS stars in a colour- 
magnitude diagram by searching for objects located above 
and to the right of the MS, since this is where one would 
expect to find very young objects. This method was first 
used by Gilmozzi et al. (1994) to identify a population of 
PMS stars outside the Milky Way, namely in the LMC cluster 
NGC 1850. More recent applications of this method include 
e.g. those of Sirianni et al. (2000), Nota et al. (2006) and 
Gouliermis et al. (2006). 

Unfortunately, this method of identification is not very reli- 
able, since the existence of an age spread and the presence in 
the same field of a considerably older population, as well as 
unaccounted patchy absorption, all tend to fill up the parame- 
ter space between the MS and the birth line, thereby thwarting 
any attempts to identify PMS stars on the basis of their effec- 
tive temperature and luminosity alone. In the specific case of 
the SN 1987 A field studied here, Panagia et al. (2000) have 
shown that the age spread is remarkable, with several gener- 
ations of young stars with ages between 1 and 150 Myr su- 
perposed on a much older field population (0.6-6Gyr). In 
such circumstances, as Figure|8] shows, broad-band photome- 
try alone would have not identified as such most of our bona- 
fide PMS stars. 

In accordance with their proximity to the MS, the major- 
ity of bona-fide PMS stars turn out to be relatively old. The 
ages of individual objects were determined by interpolating 
between the isochrones in the H-R diagram. They range from 
1 to ~ 50 Myr, with a median age of 13.5 Myr, in very good 
agreement with the estimated age of both SN 1987 A and its 
nearby Star 2 of < 13 Myr (see Scuderi et al. 1996 and ref- 
erences therein). Panagia et al. (2000) find that about 35 % 
of the stars younger than 100 Myr in this field have a typical 
age of 12 Myr, again consistent with the estimated age of the 
bona-fide PMS stars. 

The relatively old age of our PMS stars might appear at 
odds with the rather large mass accretion rates that we ob- 
tain for that population, namely ~ 10~ 8 M© yr" 1 . Accord- 
ing to current models of viscous disc evolution (see Hart- 
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FIG. 9. — The mass accretion rate as a function of stellar age for our bona- 
fide PMS stars (diamonds) is compared here with that of a similar number 
of Galactic T Tauri stars (see legend) from the work of Sicilia-Aguilar et 
al. (2006), from whom all Galactic data-points are taken (the large cross 
indicates the uncertainties as quoted in that paper). The solid line shows the 
evolution for current models of viscous disc evolution from Hartmann et al. 
(1998). Our measurements are systematically higher than the models, but 
are in agreement with the M acc values of G type stars in the Galaxy (large 
circles). 

mann et al. 1998; Muzerolle et al. 2000), the mass accre- 
tion rate of a ~ 14 Myr old PMS population should be of or- 
der 3 x 1O~ 1() M0 yr" 1 . Those models, shown in Figure|9]as a 
solid line, adequately reproduce the trend of decreasing M acc 
with stellar age observed for very low-mass Galactic T- Tauri 
stars (see e.g. Calvet et al. 2004), but under-predict our mea- 
surements (indicated by diamonds) by over an order of mag- 
nitude. 

On the other hand, recent observations of PMS stars in the 
cluster Trumpler37 by Sicilia-Aguilar et al. (2006), from 
whom the data-points in Figure|9] are taken, indicate that the 
value of M acc for — 10 -20 Myr old G-type PMS stars (large 
circles) is considerably larger, indeed of order ~ 10~ 8 M© yr, 
in excellent agreement with our results. With a median mass 
of 1.2 M Q and a median V -I colour of 0.56, our bona-fide 
PMS stars are fully consistent with the G spectral type and 
suggest that stars of higher mass have a higher M acc value 
at all ages. As mentioned in the Introduction, Muzerolle et 
al. (2003; 2005), Natta et al. (2004; 2006) and Calvet et al. 
(2004) suggest a dependence of the type M acc oc M 2 from the 
analysis of a sample of low-mass and intermediate-mass T- 
Tauri stars, albeit with a large uncertainty on the index (see 
also Clarke & Pringle 2006). A similar trend is proposed by 
Sicilia-Aguilar et al. (2006), yet with even lower statistical 
significance due to the limited mass range covered by their 
observations. The range of masses that we cover here is also 
quite limited and does not allow us to address this issue in 
detail. We will, however, do so in a forthcoming paper (De 
Marchi et al. 2008, in preparation), where we compare the ac- 
cretion rates in this field with those of a much younger and an 
order of magnitude more numerous population of PMS stars 
in the Small Magellanic Cloud. While in the SN 1987 A field 
the detection limits imposed by the WFPC2 photometry do 
not allow us to reach stars with mass below ~ 0.9 M Q , in that 
paper, based on ACS observations, we will study stars down 
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FIG. 10. — Simulated effects of a reddening increase on the mass deter- 
mination via interpolation through the evolutionary tracks. The solid lines 
are the track of Siess et al. (2000) corresponding to the masses as indicated. 
Small circles and squares show the reddening displacement, from the unred- 
dened locations (larger circles and squares), in increments of E(V — I) = 0. 1 
according to the extinction law of Scuderi et al. (1996) for this field. The 
dot-dashed line indicates the line of constant radius, while the thick dashed 
line is the same 40 Myr isochrone as in Figure[8] 
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6. SYSTEMATIC ERRORS 



The main sources of systematic uncertainty on the derived 
mass accretion rates are (/) discrepancies in the isochrones 
and evolutionary tracks, (z'z) reddening, (Hi) Ha emission gen- 
erated by sources other than the accretion process, and (z'v) the 
contribution of the nebular continuum to the colours of the 
stars. We discuss all these effects in this section. 

6.1. Models of stellar evolution 

As explained in Sections 4 and 5, the mass and age of 
PMS stars are determined by interpolation of their location 
in the H-R diagram, using model evolutionary tracks and 
isochrones as references. Apart for possible inaccuracies in 
the models input physics and errors in the interpolation, the 
largest source of systematic uncertainty on the derived mass 
and age comes from the use of models that might not prop- 
erly describe the stellar population under study (e.g. because 
of the wrong metallicity) and from differences between mod- 
els of various authors. 

For instance, if we had used Siess et al.'s (2000) tracks for 
metallicity 1/2 solar in place of those for 1/3 solar, the masses 
of our PMS objects would be systematically higher by about 
20 %. For metallicity 1/2 solar, using the tracks of Siess et 
al. (2000) or those of D'Antona & Mazzitelli (1997) results 
in masses that agree to within a few percent. On the other 
hand, when it comes to the determination of the age, using 
D'Antona & Mazzitelli's (1997) isochrones one would derive 
ages twice as young for stars colder than — 6500K, while for 
hotter stars ages would only be 20 % lower. 

In summary, systematic uncertainties of order 20 % are to 
be expected for the mass and possibly higher for the age de- 
termination. 

6.2. Reddening 



Dust extinction systematically displaces stars in the H-R 
diagram towards lower luminosities and effective tempera- 
tures. This is shown graphically in FigurjlOl for six model 
stars, three with r e g = 6 000 K (large circles) and three with 
Teft = 7000K (large squares). The effects of a reddening in- 
crease from E(V - I) = to E(V - I) = 1 on the location of 
the stars in the H-R diagram is shown by the smaller sym- 
bols, corresponding to increments of E{V-I) = 0.1. We have 
adopted here the extinction law as measured specifically in 
the field of SN 1987A by Scuderi et al. (1996). 

Since the reddening vector in the H-R plane is almost paral- 
lel to the lines of constant radius (as shown by the dot-dashed 
line for R = 10 Rq), extinction affects only marginally the es- 
timate of the stellar radius. However, at low temperatures the 
reddening vector crosses the evolutionary tracks (solid lines), 
thus leading to an underestimate of the stellar mass, if no ex- 
tinction correction is applied. In turn, this will result in a sys- 
tematically higher value of the R*/M* ratio and, therefore, in 
an overestimate of the mass accretion rate derived via Equa- 
tion|6] 

In order to better gauge the uncertainty on the value of 
R*/M* caused by unaccounted extinction, we have deter- 
mined the mass and radius of each one of the artificially red- 
dened model stars shown in Figure[l0] in exactly the same 
way in which we did that for bona-fide PMS stars, and com- 
pared those with the known input values. The result is shown 
graphically in FigureQTl where, the relation between the de- 
rived and input R^/M* ratios is shown as a function of colour 
excess E(V-I). 

Regardless of the input parameters (L, T e ff or R^/M*), it 
appears that over the range explored here the derived R*/M* 
ratio increases slowly with colour excess (see dashed line that 
best fits the median values indicated by the crosses). In partic- 
ular, underestimating the E(V-I) colour excess by — 0.4 mag 
would lead to a 30% overestimate of R*/M* (and hence of 
A^acc). This appears unlikely, considering that the best es- 
timate of the total extinction — including the Milky Way 
component — towards Star 2, one of the two bright close 
companions to SN 1987A, is E(V-I) = 0.21 (Scuderi et al. 
1996), with a la dispersion over the field of — 0.07 (Ro- 
maniello 1998; Panagia et al. 2000). According to Zaritsky 
et al. (2004), the average reddening towards cool (5500 K< 
r eft <6 500K) stars in the LMC is A v = 0.43, thus implying 
E(V-I) ~ E(B-V) ~ 0.14. Therefore, for an average LMC 
star forming region, omitting the extinction correction would 
result in a 10 % overestimate of the mass accretion rate. This 
error, albeit systematic, is smaller than most other systematic 
uncertainties and comparable to the measurement errors. 

Nonetheless, we would like to stress here that the effects 
of patchy absorption can be more severe in regions of high 
extinction. Therefore, whenever possible, one should apply 
extinction corrections to each individual stars, as it was done 
for the photometric catalogue of Romaniello et al. (2002) that 
we use here. 



6.3. Ha emission not due to the accretion process 

The underlying assumption in the determination of the mass 
accretion rate of PMS stars is that the gravitational energy L acc 
released in the accretion process goes into heating of the gas 
at the boundary layer (see e.g. Hartmann et al. 1998). In this 
case, the luminosity of the Ha emission line can provide a 
measure of the accretion energy, because it acts as a natural 



Photometric determination of the mass accretion rates of pre-main sequence stars in the SN 1987 A field 



11 



0.2 - 



0.1 



0.0 



- 


1 1 1 1 1 1 1 1 1 

• * * 

/ • 

• X * 

/ * i 


• 




s • 


X . • 
t 


X 

1 ... 1 


• 

1 . . . 1 



0.0 



0.2 0.4 
E(V-I) 



0.6 



FIG. 1 1 . — Overestimate of the R/M ratio as a function of E(V—I) when the 
extinction remains unaccounted. Th ick dots represent the R/M ratios of the 
artificially reddened stars (see Figure TlOl and the crosses their median values. 

"detector" of the luminosity released in the accretion process. 

Thus, the issue to address is whether the measured L(Ha) is 
produced only and entirely as a result of the energy released in 
the accretion process. Obviously, depending on the geometry 
and orientation of the circumstellar disc and on the specifics 
of the accretion process, some of the ionising energy could 
escape without effectively ionising the local gas. This would 
result in an underestimate of the true mass accretion rate for 
some stars, but this is true of any method that uses H emission 
lines as indicators of the accretion process and, therefore, it 
does not affect the comparison of our results with those of 
others (see Figure|9). 

One possibility that we can readily exclude is that the Ha 
emission that we detect is due to chromospheric activity. The 
total chromospheric emission of a solar-type star on the sub- 
giant branch corresponds to a luminosity of about IO^Lq 
(Ulmschneider 1979; Pasquini et al. 2000), i.e. two orders of 
magnitude less than what we measure. 

On the other hand, we must also consider the possibility that 
Ha emission might occur in discrete H knots along the line of 
sight, in places unrelated with the PMS object. Furthermore, 
hot nearby massive stars may ionise the H gas in which PMS 
stars are immersed, eventually producing Ha emission when 
the gas recombines. Both these effects would result in a over- 
estimate of M acc , but we show below that the probability of 
this happening is very small. 

As for intervening Ha emission along the line of sight, it 
would have to arise in knots of ionised H, e.g. very com- 
pact HII regions, contained within our photometric apertures. 
However, if such knots exist, they should not project only 
against stars and, therefore, narrow-band images should re- 
veal Ha emission where broad-band imaging shows no or 
very low signal (see Section 6.3). Since this is not observed, 
statistical arguments rule out this possibility as a significant 
source of the observed Ha excess emission. 

While it is possible that some of the Ha emission that we 
detect is due to diffuse nebular emission in the HII region not 
powered by the accretion process, if the emission is extended 
and uniform over an area of at least 0."5 radius around the star 
its contribution cancels out with the rest of the background 



when we perform aperture photometry (using an aperture of 
0."2 radius for the object and an annulus of 0."3-0"5 for the 
background). The subtraction would not work if the emission 
were not uniform, as for example in the case of a filament that 
projects over the star but that does not cover completely the 
background annulus. However, we have visually inspected all 
objects with excess Ha emission and excluded the few cases 
where a spurious filament of this type could be seen (see Ro- 
maniello 1998; Panagia et al. 2000). 

Another case in which the emission does not cancel out is 
when the gas ionised by an external source is associated to the 
star and completely contained within our photometric aper- 
ture. We recall here that the typical size of a circumstellar 
disc is of order 100 AU (e.g. Hartmann et al. 1998), or about 
0.02 pixel at the distance of SN 1987A, and that the median 
L(Ha) that we measure for our PMS stars is 4 x 10 31 ergs" 1 , 
or about 10~ 2 L Q (see Figure[5]l. We can then calculate how 
close a PMS stars should be to a hot, UV-bright star for the 
latter to flood the PMS circumstellar disc with enough ionis- 
ing radiation to cause at least 10 % of the Ha luminosity that 
we measure from that source. 

Let r = 100 AU be the external radius of the circumstellar 
disc and d the distance to an ionising star that illuminates the 
disc, assumed to have zero inclination (i.e. face on), so that 
4> = r 2 /(Ad 2 ) is the solid angle subtended by the disc. Let 
be the recombination coefficient to the first excited level in an 
optically thick H gas. If V is the volume of the gas and n e and 
n p are the electron and proton densities per unit volume, the 
global ionisation equilibrium equation (case B; see Baker & 
Menzel 1938) is then: 



N L 4> = (3 2 n e n p V 



(8) 



where Nl is the rate of ionising photons produced by the ex- 
ternal star. Denoting with the Ha emission coefficient for 
the same gas, the luminosity of the Ha recombination line can 
be written as: 



L(Ha) = a%%n e n p V 



(9) 



which allows us to express the luminosity of the Ha emission 
caused by the external star as a function of its ionising photon 
rate N L : 



L sx \Ha) = a 3 .2/p 2 N L cj). 



(10) 



The minimum distance beyond which the contribution of 
the ionising star to the total Ha emission becomes negligible 
(i.e. 10 % or less) is given by: 



d- = — 



r 2 a-32 10N L 



(11) 



4 p 2 UHa) 

TableQ] gives the value of 4> for a typical Ha luminosity 
of 10~ 2 Lq, as a function of the spectral type of the ionising 
star. Here we have used the compilation of Panagia (1973) 
for the appropriate stellar parameters. The table also provides 
the minimum distance d m \ n in pc and in arcsec (at the distance 
of SN 1987 A) for a circumstellar disc of radius r = 100 AU. 
We recall here that 1 arcsec corresponds to 10 WFPC2 pix- 
els. The brightest star in the field studied here is an object 
with log r e ff = 4.7 and log(L/L ) = 5.8, with an implied mass 
of ~ 60 M Q and an age less than 1 Myr (see Panagia et al. 
2000). The corresponding value of Nl, from Panagia (1973), 
is 4 x 10 49 s" 1 , thus implying d min = 0.95 pc or 38 WFPC2 
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NOTE. — The values of <j> and d mm are given for a 100 AU circumstellar 
H disc with a temperature of 10000 K. For gas temperatures of 5 000K and 
20 000 K the value of <f> is respectively higher and lower by 7.5 %, while d mm 
is 3.6% larger and smaller, respectively. 

pixel for a gas temperature of 10000K. A careful inspec- 
tion of the images shows that none of our candidate bona-fide 
PMS stars falls closer to this star than t/ m ; n . The same applies 
to the second and third brightest stars in the field, respec- 
tively with logoff = 4.7, log(L/L Q ) = 5 and logoff = 4.67, 
log(L/L ) = 4.75, as well as to all fainter objects. 

It is convenient to calculate, for each PMS object, the total 
Ha luminosity L ext (Ha) caused by all neighbouring ionising 
stars, expressed as: 

Te xt,rj , r "3,2 sr^(N L )i 

(Hd)= ^^\ x ^f (12) 

where (N£)i is the ionising photon rate of star i and di is its 
distance from the PMS object being considered. Besides the 
three brightest stars mentioned above, we have included as 
potential sources of the detected Ha emission an additional 
12 stars in the field with logoff > 4.4 and log(L/L ) > 4, as 
classified by Romaniello et al. (2002). We note that these 
estimates of L ext (Ha) are firm upper limits because we are 
assuming that (i) all disks are face-on relative to the hot stars, 
and that (ii) the projected distances on the sky are indeed the 
actual distances from the hot stars. The values of L ext (Ha) 
obtained in this way for the 133 bona-fide PMS stars span the 
range 10~ 6 - 10~ 4 Lq, or typically three orders of magnitude 
less than the median PMS Ha luminosity, as shown by the 
histogram in Figure[T2] We can, therefore, conclude that in 
this specific field the contribution of hot ionising stars to the 
measured excess Ha emission is negligible. We must stress, 
however, that the effect might be much more pronounced in 
younger star forming regions, such as the Orion Nebula or 30 
Dor. 

6.4. Nebular continuum 

Finally, consideration must be given to the role of the nebu- 
lar continuum, i.e. radiation caused by bound-free and free- 
free transitions in the gas and unrelated to the stellar photo- 
sphere. If present, nebular continuum will add to the intrin- 
sic continuum of the star, thereby affecting both the observed 
total level and the slope. This could ultimately alter the mea- 
sured broad-band colours of the source, thereby thwarting our 
attempts to infer the level of the continuum in the Ha band 
from the observed V and / magnitudes. It is important to esti- 
mate the magnitude of this effect, as our method relies on the 
intrinsic V-I colour of a star to derive its AHa (see Equa- 
tion[2]>. 

Fortunately, the contribution of the nebular continuum ap- 
pears to be insignificant. To prove this, we have assumed 
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FIG. 12. — Histogram of the ratio of observed and external Ha luminosity 
for our bona-fide PMS stars. The observed Ha luminosity is typically three to 
four orders of magnitude larger than the contribution due to external sources. 



TABLE 2 

Nebular emission and continuum of fully ionised H gas 
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NOTE. — F V ,F\ and E(Ha) are given here in erg cm 3 s Hz~' , i.e. per 
unit volume (V), unit electron density (n e ) and unit proton density (n p ). The 
luminosity can be derived by multiplying the values listed here by Vn e n p . 



a fully ionised gas of pure H, considering only bound-free 
and free-free transitions and ignoring the contribution to the 
continuum from two-photon emission (Spitzer & Greenstein 
1951). Using Osterbrock's (1989, Chapter 4) tabulations, we 
find Ha line intensity and Ha continuum fluxes of the nebular 
gas as shown in Table|2] for gas electron temperatures in the 
range 5 000 K - 20 000 K. The purely nebular W eq (Ha) ranges 
from 5 000 A to 9 000 A, or two orders of magnitude higher 
than what we measure for PMS objects. We can, therefore, 
safely conclude that the nebular contribution to the continuum 
is negligible (less than 1 %). 

Furthermore, we find that for gas temperatures in the range 
5 000 K - 1 000 K the V-I colour of the nebular continuum 
varies from 1.4 to 0.5, spanning a range typical of G-K type 
stars. Thus the effects of the nebular continuum on the V-I 
colour of PMS stars remains insignificant even for the objects 
with the highest W eq {Ha) in our sample. 

7. DISCUSSION AND CONCLUSIONS 

We have developed and successfully tested a new self- 
consistent method to reliably identify all PMS objects ac- 
tively undergoing mass accretion in a resolved stellar popula- 
tion, regardless of their age, without requiring spectroscopy. 
The method combines broad-band V and / photometry with 
narrow -band Ha imaging to: (1) identify all stars with excess 
Ha emission using as a reference template of the photospheric 
luminosity in the Ha band the average V -Ha colour of stars 
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with small photometric errors; (2) convert the excess Ha mag- 
nitude into Ha luminosity L(Ha) via the absolute photometric 
calibration of the Ha band; (3) convert, if desired, the excess 
Ha magnitude into the Ha emission equivalent width, using 
model spectra to fit the level of the continuum in the Ha band; 
(4) derive the accretion luminosity L acc from L(Ha) using a 
relationship based on recent literature values; (5) obtain the 
mass accretion rate M acc from the free-fall equation that links 
it to L acc via the stellar parameters (mass and radius), that we 
derive from the H-R diagram using pre-main sequence tracks 
for the appropriate metallicity. 

Since no spectroscopy of individual objects is needed, ob- 
servations requiring just a few hours of exposure time can 
suffice to reveal hundreds of PMS stars simultaneously, de- 
termine their Ha luminosity with an uncertainty of less than 
15 % and derive their mass accretion rates with the same ac- 
curacy as allowed by conventional spectral line analysis. 

As regards the Ha luminosity, the virtue of this new method 
is that it derives the reference luminosity of the stellar photo- 
sphere inside the specific Ha band simply from the average 
V—Ha colour of stars with the same V-I index (see Figure©. 
This is possible because, at any given time, the majority of 
the stars in a stellar field have no excess Ha emission. Ob- 
viously, the average luminosity of the stellar photosphere in 
the Ha band reflects more or less pronounced Ha absorption 
features, depending on the spectral type of the objects, and 
precisely for this reason it defines the reference luminosity 
with respect to which any excess emission should be sought 
and measured. The uncertainty associated with the derived 
L(Ha) is dominated by statistical errors in the photometry and 
does not exceed 15 % for our 4a selection (see Section 3.1). 
Note that this approach is also applicable to populations com- 
prising mostly very young objects, such as the Orion Nebula, 
since excess Ha emission is found in only about one third 
of the PMS stars at any given time and, therefore, the major- 
ity of the very young objects will not have a conspicuous Ha 
emission, thus correctly defining the reference template (see 
Section 3.1). 

In spite of the good accuracy on L(Ha), deriving the mass 
accretion rate requires the use of a rather uncertain relation- 
ship linking L(Ha) to the accretion luminosity (see Section 4 
and also Dahm 2008). The uncertainty on the relationship can 
have a systematic effect of up to a factor of 3 on the value of 
M acc . This is the unfortunate consequence of our incomplete 
understanding of the physics of the accretion process, plagued 
by the lack of direct measurements of the bolometric accre- 
tion luminosity L acc , as well as by the paucity of simultaneous 
measurements of the hot continuum excess emission and of 
the emission line fluxes in Ha or other lines (e.g. Pa/3 and 
Br7). Although annoying, the systematic uncertainty of up to 
a factor of 3 on M acc is similar to that plaguing age determi- 
nations based on the comparison with theoretical isochrones, 
which can reach a factor of 2 (see Section 6.1). 

On the other hand, these systematics affect in the same way 
all determinations of M acc that make use of emission lines di- 
agnostics, including those based on a detailed study of the 
profile and intensity of the Ha line (e.g. Muzerolle et al. 
1998b), since their calibration relies on the same measure- 
ments. For this reason, the relative values of M acc that we 
obtain are trustworthy and the comparison of our results with 
those obtained with other methods remains meaningful. 

As discussed in Sections 3 and 4, as an application of our 
method, we have studied the field of SN 1987 A, where we 
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FIG. 13. — Fraction of PMS objects with respect to all stars, as a function 
of the distance (in units of WFPC2 pixel or Of 1) from the baricentre of the 
ionising stars. 

identify 133 bona-fide PMS stars, based on their large Ha ex- 
cess. The median value of the mass accretion rate that we 
find for these stars, namely 2.9 x 10~ 8 M Q yr" 1 , is in very 
good agreement with the median value of 2.4 x 10~ 8 M Q yr" 1 
found by Romaniello et al. (2004) from the analysis of the 
distribution of the (U-B) colour excess in the same field. 
While we have identified 133 bona-fide PMS stars with an 
Ha excess above the 4 a level, Romaniello et al. conclude 
that 765 stars with an excess U -band luminosity larger than 
~ 0.035 L© could be bona-fide PMS candidates, albeit at the 
1 .5 a level. 

The approach followed by Romaniello et al. (2004) is sta- 
tistical in nature and, while providing a solid estimate of the 
average mass accretion rate, it does not allow one to reliably 
identify individual PMS stars. Therefore, Romaniello et al. 
turned to the Ha excess emission (namely R- Ha) as a more 
direct indicator of ongoing accretion and found a strong cor- 
relation between the f/-band continuum excess and W eq (Ha) 
for stars with U-B > 0.5 and W ecj (Ha) > 3 A. Not surpris- 
ingly, the more stringent threshold on the equivalent width 
adopted in the present study, i.e. W eq (Ha) > 10 A, results in a 
smaller number of identified of bona-fide PMS stars, but also 
in a much smaller uncertainty. Indeed, the method that we 
present here is more direct than the one based on the f/-band 
excess, since it does not require the use of model atmospheres 
to predict the level of the intrinsic photospheric continuum in 
the U band from the observed optical colours, and as such it 
results in a more reliable determination of the accretion lumi- 
nosity. 

Our method, therefore, makes it possible for the first time 
to carry out a quantitative and systematic study of the proper- 
ties of a large number of PMS stars, including those in dense 
and/or distant star forming regions. This allows us to tackle 
some still open questions as to the evolution of both the stars 
and their discs. For instance, in the field of SN 1987A, we 
have studied how the mass accretion rate changes with time 
as our bona-fide PMS stars approach their main sequence. We 
find that M acc decreases more slowly with time than what 
is observed for low-mass T-Tauri stars in the Galaxy (see 
Section 5) and predicted by models of viscous disc evolution 
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FIG. 14. — Map of the distribution of all bona-fide PMS objects with re- 
spect to the ionising stars (blue star symbols), whose baricentre is indicated 
by the large yellow square. Red squares correspond to stars with L(Ha) > 
2 X 10~ 2 Lq, white squares indicate objects with L(Ha) < 8 X 10~ 3 Lq, and 
green diamonds are used for intermediate values. Note the paucity of PMS 
objects with high Ha luminosity near the baricentre of the ionising stars. 



(Hartmann et al. 1998). However, our results are in excellent 
agreement with the (so far limited) measurements of the mass 
accretion rate of G type stars in the Milky Way (e.g. Sicilia- 
Aguilar et al. 2006), thus lending credit to the hypothesis that 
Mace might depend on the stellar mass. 

Moreover, our method allows us to investigate whether and 
how the accretion process is affected by the chemical compo- 
sition and density of the molecular clouds or by the proximity 
of hot, massive ionising stars. We will address the effects of 
metallicity in a forthcoming paper (De Marchi et al., in prepa- 
ration) devoted to a comparison of star formation in the Small 
and Large Magellanic Clouds. As for the effect of hot massive 
stars, we have studied the distribution of the 133 bona-fide 
PMS objects with respect to the baricentre of the 15 hottest 
stars in the field, described in Section 6.3. We find a clear 
anti-correlation, shown in Figure[13] between the frequency 
of the identified PMS stars and their distance from the ionis- 
ing objects (marked with blue star symbols in Figure[l4l. This 
effect is not seen for non-PMS objects of similar brightness, 
whose distribution is uniform over the whole field, confirming 
that the trend is not due to problems in detecting faint objects 
near the brightest stars. The trend remains the same whether 
the baricentre of the hot stars is computed using as weight 
their bolometric luminosity or their ionising flux. This anti- 
correlation is quite remarkable in several ways, as we discuss 
below. 

In Figure[l4]we show a schematic map of the distribution of 
all 133 PMS objects with respect to the 15 ionising stars. The 
baricentre of the ionising flux is shown by the large square 
and concentric annuli are drawn from it, in increments of 35". 
For reference, a red ellipse marks the position of SN 1987 A. 
Symbols of different colours indicate PMS stars with differ- 
ent Ha luminosity, from red for L(Ha) > 2 x 10~ 2 L Q to white 
for < 8 x 10~ 3 L Q , with green diamonds corresponding to in- 



termediate values. An inspection to Figure[l4] immediately 
reveals two important facts: (?) there are very few PMS stars 
near the baricentre of the massive objects and (ii) their Ha lu- 
minosity is systematically lower than that of PMS stars farther 
away. 

Panagia et al. (2000) had already pointed out that there is no 
overdensity of young PMS objects around hot massive stars 
in this field. This shows that low-mass and massive stars do 
not follow each other, suggesting that their formation mecha- 
nisms are separate and distinct. But what Figure[14] addition- 
ally shows is that there are very few PMS stars with strong Ha 
luminosity (red squares) around the young massive stars, with 
only 20 % of them within a 70 arcsec radius of the baricentre. 
Conversely, stars with low Ha luminosity (white squares) are 
distributed rather uniformly across the field. 

With a median age of 13.5 Myr, most PMS objects are much 
older than the 1-2 Myr old hot stars and must belong to a 
previous generation. Thus, no spatial relationship (whether a 
correlation or anti-correlation) should be expected between 
the distribution of the two types of objects. On the other 
hand, while we cannot exclude that originally fewer PMS 
stars formed in the region now populated by young massive 
stars, there are indications that the local paucity of PMS ob- 
jects is only apparent and caused by the presence of the hot 
young stars. 

The fact that PMS objects near the youngest hot stars are 
both less numerous and fainter in Ha emission suggests that 
their circumstellar discs have been considerably eroded and 
made less efficient by enhanced photo-evaporation. We have 
shown earlier (see Section 6.3 and Figure[T2b that the instan- 
taneous contribution to the Ha luminosity from the ionising 
radiation of nearby massive stars is negligible. The results 
shown in Figure[l4] confirms that this is the case. If it were 
not and the observed Ha luminosity were mostly due to fluo- 
rescence of the circumstellar discs under the effect of the ion- 
ising radiation of massive stars, the radial distribution should 
be opposite to the one observed, with brighter Ha luminosi- 
ties for PMS stars closer to the baricentre of the youngest hot 
stars. Thus, if the apparent paucity of PMS objects in their 
vicinity is the result of photo-evaporation, the latter must be 
mostly due to the effects of non-ionising far ultraviolet (FUV) 
radiation. When integrated over the lifetime of massive stars 
(~ 2 Myr), their FUV radiation can be quite significant for 
the disruption of nearby circumstellar discs. In particular, this 
also means that there might be many more objects still in their 
PMS phase, close to the hot stars, but we have no way to iden- 
tify them as such, since the photo-evaporation of their discs 
has stopped the accretion process. 

Recent theoretical studies of the erosion of circumstellar 
discs via photo-evaporation due to external FUV radiation 
(e.g. Adams et al. 2004) or to the PMS star itself (e.g. Alexan- 
der, Clarke & Pringle 2006a; 2006b) suggest timescales of 
order 8-10 Myr for the complete dissolution of the discs. Al- 
though this is the right order of magnitude for the timescale of 
the process that we see, there are still some discrepancies. For 
instance, it is clear that if discs were to completely dissolve 
under the effect of their own central star within 5-10 Myr, 
as suggested by Alexander et al. (2006b), the majority of our 
PMS objects, with a median age of 13.5 Myr, should not show 
any sign of accretion. Similarly, with a typical age of 2 Myr, 
the hottest stars in our field would have not had time to de- 
posit enough UV photons on to the PMS discs to seriously 
erode them, if the process is to take of order 10 Myr (Adams 
et al. 2004). While recently Clarke (2007) suggested that the 
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effect of nearby massive young stars may be stronger, leading 
to disc exhaustion timescales of order 2 Myr, her calculations 
apply to the immediate vicinity (< 1 pc) of very massive stars, 
such as the 06 star 0\C in Orion. The density of massive stars 
in the field of SN 1987 A is much lower than that present in the 
Orion Nebula and the distances to our PMS objects typically 
an order of magnitude larger than those considered by Clarke 
(2007). 

Nevertheless, while it is clear that the theory of disc photo- 
evaporation does not yet capture all the details, we believe that 
it can greatly benefit from the observations presented here. In 
particular, it will be necessary to explain which mechanisms 
allow isolated circumstellar discs to feed their central PMS 
stars for well over 20 Myr, at least in a low-metallicity, low- 
density environment such as the field of the LMC. At the same 
time, the models should justify the quick (< 3 Myr) demise 
of the PMS discs within a (projected) radius of ~ lOpc of a 
handful of late O stars. Extending our investigation, as we 
plan to do, to more and denser regions, such as NGC 1850 
and 30 Dor in the Large Magellanic Cloud and NGC 346 and 
NGC 602 in the Small Magellanic Cloud, will offer increased 
statistics and a much wider variety of star formation environ- 
ments. 
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the WFPC2 filter system based on all the observed spectra in the BPGS atlas, 
with solar metallicity. The range over which the least square fit is carried out 
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APPENDIX 

Over the years, a large number of observations of stel- 
lar populations have been collected with the HST/WFPC2 
and HST/ACS through the V-like filters F555W and F606W 
and the /-like filter F814W. We have therefore computed 
colour relationships between the magnitude in the R-like band 
F675W (niffis) and the V-7-like colours mm -mm and 
'«606 _ 'W8i4- To this aim, we have used both observed and the- 
oretical stellar spectra, taken respectively from the Bruzual- 
Persson-Gunn-Stryker (BPGS) Spectrophotometry Atlas (see 
Gunn & Stryker 1983) and from the stellar atmosphere mod- 
els of Bessell et al. (1998). 

The BPGS atlas contains nearby bright dwarfs and giants 
of near-solar metallicity, whereas from the library of Bessell 
et al. (1998) we have selected specifically those models with 
metallicity index [M/H] = -0.5, as appropriate for the LMC 
(Dufour 1984), surface gravity g = 4.5 and effective tempera- 
tures in the range 3 500 K < T eii < 40 000 K. Using the stan- 
dard HST synthetic photometry package Synphot (Laidler et 
al. 2008), we have folded these spectra through the instrumen- 
tal response of the WFPC2, the ACS/WFC, the ACS/HRC 
and the WFC3/UVIS cameras, in order to obtain the ex- 
pected magnitudes in the bands listed above. For comparison 
with ground-based photometry, we have also calculated the 
expected magnitudes through the Johnsn-Cousin V,R and / 
bands. 

We show in Figure[l5] for the WFPC2 filter system, the 
predicted 7W555-7W675 colour as a function of 7W555-7W814 for 
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FIG. 16. — Relationship between the 111555— "7675 and "7555 — "1814 colours 
of the WFPC2 filter system based on model spectra of Bessell et al. (1998) 
for metallicity [M/H] = —0.5. The range over which the least square fit is 
carried out (4 000 K < r cff < 10 000 K) is indicated by filled circles. 

all the stars in the BPGS atlas. The distribution is remark- 
ably linear and the dashed line represents the best fit in the 
colour range < 771555-771814 < 1.5 (filled symbols), corre- 
sponding to the relationship 771575 = 0.34 x 771555+O.66 x 7«si4. 
The colour-colour diagram in the same bands based on the 
models of Bessell et al. (1998) is shown in Figure[T6land the 
relationship is practically indistinguishable, namely 777675 = 
0.35 x 777555 + O.65 x 7«8i4. In this case, the dashed line shows 
the least-square fit to stars in the range 4000K < T e ff < 
10000K, marked by the filled symbols. As already noted in 
Section 2, it appears that neither the luminosity class nor the 
metallicity has a noticeable impact on the colour relationship 
in these bands. 

This is also the case when we consider the tt?606 _ "isi4 
colour, for which the relationship with 7W675 is given in Ta- 
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TABLE 3 

Predicted colour relations between R or Hcf and V- 



I FOR VARIOUS PHOTOMETRIC SYSTEMS. 



System 


Band 


BPGS, [M/H] = 


Bessell etal., [Af/ff] =- 


0.5 


Bessell et al., [M/H] = - 


-1 


WFPC2 


'"675 
m 61s 

m 656 
<K 


0.342 x m 5 5 5 + 0.65 8 x m 8 i 4 
0.510 x m 606 + 0.490 x m 814 


0.354 x m55j +0.646 x m 8 i 4 
0.510 x m 60 6 +0.490 x m 814 
0.409 x fn 555 + 0.591 x m 814 
0.593 x m 606 + 0.407 X m 814 


-0.351 
-0.353 


0.353 x »J555 + 0.647 x m 8 i 4 
0.508 x m 606 + 0.492 x m 814 
0.411 XM555 +0.589 xm 8M 
0.595 X m 60 6 + 0.405 X m 814 


-0.350 
-0.351 


ACS/ 
WFC 


»<625 
m 625 
<5S 
<K 


0.496 X ms55 +0.504 X m 8 i 4 
0.708 x m 6 06 + 0.292 x m 8 i 4 


0.499 X OT555 +0.501 X m 814 
0.703 x m 606 + 0.297 x m 814 
0.381 x m 55 5 + 0.619 x m 8M 
0.546 x m 606 + 0.454 x m 814 


-0.156 
-0.160 


0.503 X JM555 +0.497 X m 814 
0.708 xm 60 6 + 0.292 x m 814 
0.386 x msss +0.614 x m 8 i 4 
0.552 X m 60b + 0.448 X m 814 


-0.158 
-0.164 


ACS/ 
HRC 


"«625 
™625 
m M8 
< S 


0.506 x m 5 55 +0.494 x m 8 i 4 
0.707 xm 606 + 0.293 x m 814 


0.502 x OT555 +0.498 x m 814 
0.697 x m 606 + 0.303 x m 814 
0.385 x m555 + 0.615 x m 8 i 4 
0.542 x m 60 6 + 0.45 8 X m 814 


-0.160 
-0.166 


0.511 X m 55 5 + 0.489 X m 814 
0.704 x m 60 6 + 0.296 x m 814 
0.390 x ms55 +0.610 x m 8 i 4 
0.544 x m 60b + 0.456 x m 814 


-0.160 
-0.165 


WFC3/ 
UVIS 


'"625 
"1625 
m 656 
<K 


0.504 X m 55 5 +0.496 X m 8 i 4 
0.742 X m 60 6 + 0.25 8 X m 8w 


0.510 X m 55 5 + 0.490 X m 814 
0.734 x m 606 + 0.266 x m 814 
0.387 x m 5 5 5 + 0.613 x m 814 
0.558 x m 60 6 + 0.442 X m 814 


-0.405 
-0.406 


0.517 X m 5 5 5 + 0.48 3 X m 814 
0.741 x m 606 + 0.259 x m 8M 
0.389 xm 55 5 + 0.611 x m 814 
0.560 X H7606 + 0.440 X m 814 


-0.408 
-0.409 


Johnson- 


R 


0.451 X V + 0.549x7 


0.465 X V +0.535 xl 




0.473 xV + 0.527 x/ 




Cousin 














Note. 


— The typical 1 a 


uncertainty on the slope of the regression 


fit is 0.003. 









1.5 



1.0 



0.5 - 



0.0 



m c =0.409 m +0.591 m -0.351 ^° 

656 555 814 , 



a' 



9' 



-0.5 



0.0 



0.5 
m 



1.0 



1.5 



2.0 




FIG. 17. — Relationship between the 'M555 - ffjjj 56 and m^ — m 8 i 4 colours 
of the WFPC2 filter system based on model spectra of Bessell et al. (1998), 
for metallicity [M/H] = — 0.5, in which all spectral features in the range 6 500 
- 6 620 A have been replaced by the continuum. The range over which the 
least square fit is carried out (4 000 K < T efl < 10 000 K) is indicated by filled 
circles. 

ble|3] In the same table we also list the results for the 
ACSAVFC, ACS/HRC and WFC3/UVIS cameras, as well as 
those for the standard Johnson-Cousin system normally used 
from the ground (note that the R and / bands in this system do 
not coincide with the R and / bands of Johnson's original sys- 
tem, as mentioned in Section 2). As expected, the coefficients 
in the colour relationship vary considerably with the band, and 
differences exist also between the same band on various HST 
cameras (or, at least, bands with the same name). 

As briefly discussed in Section 3, the magnitude corre- 
sponding to the level of the continuum in the Ha band can be 
derived by means of properly validated model atmospheres, 
such as those of Bessell et al. (1998) that we have shown to 
reproduce quite accurately the observed broad-band colours 
(see Figure|2]i. This means that, even though the models might 
lack the spectral resolution necessary to realistically repro- 



FlG. 18. — Comparison between the Ha absorption equivalent widths in the 
catalogue of Ducati (1981; crosses indicate dwarfs, squares are for giants) and 
those derived by applying Equationf4]to the spectral models of Munari et al. 
(2005) for solar metallicity (solid line). 

duce spectral lines (hence the small discrepancy between the 
dashed lines and the squares in Figure|2]i, the level of the con- 
tinuum can be trusted. We have, therefore, fitted the con- 
tinuum in the Bessell et al. (1998) spectra over the range 
6 500 -6 620 A and have replaced the actual spectrum in that 
range with the interpolated continuum. Finally, we have used 
Synphot to fold the model spectra, modified in this way, 
through the instrumental set-up and derive in this way the 
continuum Ha magnitudes (Ha c ) for the various Ha filters 
on-board the HST, as a function of the spectral type or T e s of 
the star. 

With this information it is now possible to build colour- 
colour diagrams like those of Figure[l5] and Q~6l and look for 
a colour relationship between Ha c and V -I. In FigurefTTI 
we show the specific case of the WFPC2. The OT555 - m c 656 
colour varies very smoothly and very slowly with the effec- 
tive temperature and is a remarkably linear function of the 
?W555 - m8i4 colour over the whole range, and particularly for 
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TABLE 4 

Rectangular width and reference wavelength of past and 
present ha filters installed on board the hst. 



Instrument 


Filter 


RW [A] 


Aref 


WFPC2 


F656N 


28.35 


6563.8 


ACS/WFC 


F658N 


74.96 


6584.3 


ACS/HRC 


F658N 


74.96 


6583.9 


WFC3 


F656N 


17.68 


6561.4 




F658N 


27.61 


6585.0 



stars with 4 000 K < r efl < 10000 K (filled circles) with which 
we are mostly concerned in this work. Similar relationships 
exist for other bands and photometric systems and they are 
given in Table[3] We note that the constant 0.351 in the colour 
relationship accounts for the conspicuous Ha absorption line 
in the spectrum of Vega, which defines the standard for the 
VEGAMAG magnitude system. 

As indicated in Section 3, the difference between the ob- 
served Ha magnitude and Ha c can be converted to the equiv- 
alent width E eq (Ha) by means of Equation[4] To prove this, 
we have applied Equation|4]to the model atmospheres of Mu- 
nari et al. (2005), which provide spectra with a resolving 
power of 1 1 500 and a 1 A sampling, and compared the re- 
sults to literature measurements of Ha equivalent widths (in 
absorption). 
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